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In a previous paper we developed a form of variational trial wave function and applied it to 
van der Waals clusters: five or less atoms of Ar and Ne modeled by the Lennard- Jones potential. In 
addition, we tested the trial functions for a hypothetical, light atom resembling Ne but with only 
half its mass. We did not study atoms such as He* with larger de Boer parameters, i.e., systems in 
which the zero point energy plays a more important role relative to the potential energy. This is the 
main purpose of the present paper. In fact, we study clusters to the very limit where the zero-point 
energy destroys the ground state as a bound state. A simple picture of this un-binding transition 
predicts the power law with which the energy vanishes as the de Boer parameter approaches its 
critical value and the power of the divergence of the the size of the clusters in this limit. Our 
numerical results are in agreement with these predictions. 



I. INTRODUCTION 



We consider clusters of bosonic Lennard-Jones atoms for which we present two sets of results. Firstly, there are 
improved estimates of the ground state energies of systems studied previously with variational Monte Carlo Q. The 
improvements were obtained with a modified diffusion Monte Carlo algorithm similar to Ref. j^. Secondly, 
we study the behavior of the clusters for small masses. That is, in reduced units such that the Lennard-Jones pair 
potential has the form r"^^ — 2r^^, the only independent parameter in the Schrodinger equation is the reduced inverse 
mass 'm~^ , a quantity proportional to the square of the de Boer parameter. As the de Boer parameter increases, the 
relative importance of the zero-point energy and the ground state energy of a cluster increase, as does its size. At a 
critical value of the de Boer parameter the ground state energy Eq vanishes while average cluster size (r) , as defined 
below, diverges and the cluster ceases to exist in a bound state. 

For the simple case of a dimer one can show that 

Eo Iminio ^ Am^, 

{r) Uimo ~ Am"\ (1) 

where Am ~ m — rUc with rric the critical value of the reduced mass. Note that rUc depends on TV, the number of 
atoms in the cluster and is expected to be a monotically decreasing function of N. The mathematical mechanism that 
yields Eqs. (|l]) is the following. Two scattering states forming a complex conjugate pair merge at zero momentum to 
produce two states with "complex momentum" : a physically acceptable bound state and a state with unacceptable 
behavior at infinity. This mechanism is probably not limited to the dimer and it is quite plausible that Eqs. (|^) apply 
in general to clusters of any finite size. 



II. RESULTS 

Table || shows the comparison between the ground state energy estimates obtained by using variational Monte 
Carlo |jl| and our improved diffusion Monte Carlo algorithm. Results obtained by variational Monte Carlo suffer 
from a systematic bias, i.e., if we denote by Et = {■iPt\'H\4>t) , the variational estimate obtained a given normalized 
trial state ItJjt), one has Eq < -Et, where the equality holds only if the trial function is the exact ground state wave 
function. 

If one defines 

= {^T\in ^ Eo fl^Pr) (2) 
the following inequality holds (see Ref. ^ for details and references): 

0<Et-Eo<^-^, (3) 
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TABLE I. Estimates of the ground state energies Eq for noble gases Ar and Ne and hypothetical lighter particle ^-Ne obtained by using 
improved diffusion algorithm compared with estimates St taken from Ref. [Q. Standard errors in the last digit are given in parentheses. 
Estimates of the relative errors, as described in text, are given by Q' and Q" . Missing values indicate cases where the statistical errors 
are smaller than the errors due to numerical differentiation. 
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where Ei is the energy of the first, totally symmetric excited state. To estimate the number of correct digits in the 
variational estimate of the ground state and to ascertain how good a bound inequality (|^) is we introduce the following 
quantities: 
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Q" = 



log 
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{El — Et){Et — Eq) 



The results are shown in Table |[ Quite remarkably, the bound given in Eq. (^) is very tight. 

Next we compare the behavior of the ground state energy on mass with the behavior predicted by Eqs. (|l]). Fig. |l| 
shows the energy as a function of the de Boer parameter for clusters of sizes iV = 3, 4 and 5. The energy has 
been normalized by dividing by the classical ground state energy and we note that the linear behavior for small de 
Boer parameter, i.e., large mass, follows from the harmonic approximation. Fig. ^ explicitly shows the data in a 
double-logarithmic plot of the ground state energy vs the deviations from the respective critical points. 

The average size of the clusters was expressed in terms of r.y, the distance between atoms i and j, using the 
following two definitions: (1) r ^ 2{r,j)/N{N - 1) (2) R = {2{rfj)/N{N - 1))^/^. Double logarithmic plots of the 
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FIG. 2. Log-log plot of the normalized ground state energy -Bo/^'clas f^- the deviation from the un-binding transition, m — rrtc. 



average size vs the deviations from the respective critical points are shown in Figs. ^ and ^ It should be noted that 
in contrast to the diffusion Monte Carlo estimates of the energy, the estimates of the average size are biased in the 
sense that for any operator A that does not commute with the hamiltonian diffusion Monte Carlo yields the matrix 
element (V'tI^IV'o) rather than the ground state expectation value (V'ol^lV'o)- The values shown in the figures were 
obtained by linear extrapolation from (i/)t|^|V't) and (i/'tI^IV'o)- Irregularities in the quality of the trial functions 
are presumably responsible for the corresponding irregularities in the average cluster sizes for the smallest masses in 
the N = 5 case. 



III. CONCLUSIONS 



We presented estimates of the ground state energies of van der Waals clusters in a wide range of masses. The results 
show that wave functions introduced in Ref. provide good trial function in the whole range from the classical limit 
to the ultra-quantum limit, where the clusters cease to form a bound state. Results for the behavior in the vicinity of 
the un-binding transition corroborates Eqs. (eq. critical). The agreement for the case of the energy is very good; for 
the divergence in the size dependence the average cluster size of the clusters closest to the transition was not quite 
large enough to show the asymptotic "critical behavior." 
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FIG. 3. Log- log plot of the average cluster size r as defined in the text vs. the deviation from the un-binding transition, m — rric- 
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FIG. 4. Log- log plot of the gyration radius R vs. the deviation from the un-binding transition, m — rric. 
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